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Viscosity Sum Rules at Large Scattering Lengths 

Walter D. Goldberger and Zuhair U. Khandker^ 
^ Department of Physics, Yale University, New Haven, CT 06520 

We use the operator product expansion (OPE) and dispersion relations to obtain new model- 
independent "Borel-resummed" sum rules for both shear and bulk viscosity of many-body systems 
of spin-1/2 fermions with predominantly short range S-wave interactions. These sum rules relate 
Gaussian weights of the frequency- dependent viscosities to the Tan contact parameter C{a). Our 
results are valid for arbitrary values of the scattering length a, but receive small corrections from 
operators of dimension A > 5 in the OPE, and can be used to study transport properties in the 
vicinity of the a — ^ co fixed point. In particular, we find that the exact dependence of the shear 
viscosity sum rule on scattering length is controlled by the function C{a). The sum rules that 
we obtain depend on a frequency scale ujq that can be optimized to maximize their overlap with 
low-energy data. 



INTRODUCTION 



^ \ Fermions with tuned short-range interactions, resulting in large S'-wave scattering length, seem to describe a wide 

j^JT)' range of physical phenomena arising in diverse areas of physics, from nuclear and hadronic processes to table top 
I ' experiments in condensed matter and atomic physics. This is because the unitary limit, in which the scattering length 
(~| . a diverges, is governed by an ultraviolet fixed point of the renormalization group in which the dynamics becomes 
d ' universal. Progress over the last few years in manipulating samples of cold atoms tuned to a Feshbach resonance has 
^ „ made it possible to study the many-body dynamics of these systems in a detailed and controlled manner. 

One of the ways in which universality manifests itself in these systems was uncovered in the work of S. Tan [lH3|, 
which derives a set of relations that hold true when evaluated in any quantum state, and which are fixed, for a — >■ oo, 
in terms of few-body dynamics at short length scales. It was subsequently pointed out by Braaten and Platter [j| 
^ that the Tan relations can be derived from the operator product expansion (OPE) Q applied to the quantum field 

i-rt theory realization of fermions near unitarity, 



o : . _ Y- .,,t / -^ , V' 



£ = ^ 7/^t z9t + — ) Va + g44i'ii^2. (1) 



Indeed, the existence of such universal relations in field theory is natural from the point of view of the existence 
^^ , of an ultraviolet fixedpoint a ^ oo with enhanced non-relativistic conformal (Schrodinger) symmetry mWi- Since 
J:^ the original work of [ijs, a large class of universal relations have been obtained using OPE methods |10l - [l3l | or 
otherwise [14| - [l9| . See 20| for a recent review of these developments. 

In this paper, we report on a new set of sum rules that relate weighted integrals of the frequency-dependent shear 

and bulk viscosities Rery(u;), ReC(a;) of a many-body system near unitarity to the contact parameter C introduced 

in [l|. In order to obtain these sum rules, we first use the OPE to fix the w -> oo asymptotic behavior of the Feynman 

(time ordered) two-point correlator of the traceless part of the stress tensor of the theory defined by Eq. ^ . (The 

T— I , analogous OPE for the particle number density is known [l^] , and that of the particle number current is presented 

L| in the appendix as an independent check of results.) By standard Kubo formulae, together with analyticity, these 

. ^H , OPEs give the tail behavior of the viscosities at w — >■ oo. The result for the tail of Re r]{uj) agrees with the previously 

j^ ■ obtained expression in ref. [19], while the asymptotic bulk viscosity is new (see also refs. [lE: [2l| for previous work 

?H ] on the viscosity of unitary gases). Unlike ref. [13], however, our expressions are valid for arbitrary values of the 

parameter a^/muj. They are accurate up to calculable corrections to the OPE from operators of dimension A > 5 

(e.g., three-body contact operators as well as finite range effects). Finally, in order to obtain results that can be 

compared to low-energy experiments, we use the methods introduced in ref. [12| (themselves based on the QCD sum 

rule technology of [22]), which employ a Borel resummation of the OPE to derive sum rules that can be made to 

overlap with the data by adjusting a free parameter ujq. 

This paper is organized as follows. Sec. |lT] collects several previous results that will be used in our derivation of 
sum rules, including the relevant Kubo formulae in sec. Ill Al and a review of the interplay between dispersion relations 
and the OPE in sec. IIIBI In sec. IIII Al we compute the leading term for the OPE of the traceless-part of the stress 
tensor with itself, working to all orders in the coupling constant g of Eg. ([IT) but to leading order (dimension A = 4) 
in the operator expansion. This, together with the previous results of [lj| on the structure factor S{u},q) allows us 
to obtain in sec. IIII B1 the asymptotic properties of both bulk and shear viscosity as well as sum rules to all orders in 
a. The scattering length dependence of the shear viscosity results is given exactly by that of the contact C{a), while 



the expressions involving bulk viscosity vanish as a~^ in the limit of infinite scattering length. In appendix]^ we give 
the full OPE for the longitudinal and transverse parts of the current two-point function, which gives an independent 
derivation of our results and serves as a consistency check. We conclude in sec. IIVI 

II. PRELIMINARIES 

A. Hydrodynamic Relations 

The connection between hydrodynamics and the microscopic theory of Eq. ([T]) is given by standard Kubo formulae. 
We will focus on viscosity, which at the macroscopic level is defined in terms of a small gradient expansion of the 
stress tensor T^^{x) : 

Tij{x) = P{x)Sij + p{x)ui{x)uj{x) - (T^j{x), (2) 

where x — {t,x). In this equation P{x) is the pressure, p{x) is mass density, and u{x) is the fluid velocity. The 
dissipative part of the stress tensor, in the large wavelength limit, is given by 



(Jij = V 



diUj + djUi - -Sij{\/ ■ u) 



+ C<5„-(V •«) + •••, (3) 



where d is the number of spatial dimensions and • • • denotes higher-order gradients of u. Thus the long wavelength 
fluid transport properties are encoded in terms of two transport coefRcients, namely the shear viscosity rj and bulk 
viscosity C. 

At the microscopic level, 77 and C are the zero-frequency limits of particular linear response functions associated 
with the quantum theory. Specifically, for a non-relativistic theory with conserved particle number J^ = (n^x), J{x)), 
consider the response function 

x'j,'{q) = iJd''+'xe^'^-^0{t){[r{x),J^{Q)]), (4) 

with the notation x^ = (f , x), q^ — (w, q), q ■ x = wt — q ■ x, and d'^'^^x = dtd'^x. Decomposing this Green's function 
into longitudinal and transverse components 



X^iq) - XLiqf-^+XT{q) (S^' - ^ ) , (5) 



q \ q 

one defines [2j| frequency dependent viscosities 77(0;) and C{uj) by 

t2 



77(0;) = hm — XT(w,q), 6) 



2d -2 



C(cj) H — 77(w) = Imi — - — XL{^,q)- (7) 



With these conventions, the dissipative response is encoded in the real parts of the complex functions 77(0;) and C('^)i 



2 

Re77(a;) = lim — ^-Imxrlw, g), (8) 

q-^o q'^ 



2d~2 m^w, 

— — Ke7;(a;j — hm —^^ 

d g-^O (7^ 



ReC(w)H -^ — Re?7(a;) = \m)^^^lv[\XL{^,q)- (9) 



The constants 77 and C, are the zero-frequency limits of these functions: 77 = Re 77(0; = 0), C = ReCC^^ — 0). 

In addition to the number density and current (77, J{x)), hydrodynamic transport is characterized by several other 



conserved quantities, namely energy density T'^'^ = H, energy flux T**' and stress T*^. These obey continuity equations 

dtn + da'^0, (10) 

dtT°° + d,T'° ^ 0, (11) 

mdtr + djT^' ^ 0. (12) 

As a result Re ri{io) and Re C{^) can also be obtained from the stress tensor correlator 

x)^'"(g) = »/d'^+'xe^«-0(i)([r^"(x),r^-^(O)]). (13) 

Thus we can write, 

Re?7(a;) = lim — Inix^^'^^(w,q), (14) 

ReC(i^) H J — Re?7(a;) = lim —lrRx^'^^{i^,q)- (15) 

a (f— >0 uj 

Additionally, conservation of number current relates the longitudinal function xl(w, g) to the correlator of the number 
density "'^ 

XL{u^,q) ^'-^ J d''+'xe^'^-^e{t){[n{x),nm), (16) 

and therefore 

9(7—9 Trm'^in'^ 

ReC(a;) + f^Rer,{u;) = hm ^^5(c.,g), (17) 

where S{u!,q) is the dynamic structure factor. 

In this paper we will find it most convenient to use Eqns. (J14p . (1171) to obtain viscosity sum rules. In order to obtain 
hydrodynamic transport coefficients for the theory defined by Eq. ([T]) we need operator expressions for the conserved 
currents. The particle number density and current operators are 

n(x) = il}H{x), (18) 

J{x) = ~-^^t ^ ^. (19) 

Am 

For the purposes of this paper, it is sufficient to define the stress-tensor as the spatial part of the Noether energy- 
momentum tensor T^^^ associated with invariance under time plus space translations, whose components are [25| 
energy density, 

H = TOo = ^|VV|'-gVlV'IV'iV'2, (20) 



energy flux. 



and stress tensor 



r^" = - ^ [{dt^^) (a, V) + (5.^^) {dt^)] , (21) 



T'' = ^ [{^^i'^) idj^) + (a,^0 (a^^)] „ s'j 



— V2(V'V^)+3V'IV'|V'iV'2 



(22) 
The remaining component is the momentum density, which obeys T'^* = rnJ' as required by Galilean invariance. 



^ This relation holds modulo Schwinger terms which arc analytic in momentum space and therefore do not contribute to dissipative 
response. 



B. OPE and Sum Rules 

The unitary limit is characterized by strong coupHng, and thus the problem of computing the functions ri{uj) and 
C(uj) in the many-body ground state is generally not analytically tractable. However, as discussed in refs. 0, Hol - 
[l2| . it is possible to obtain some analytic results in the w — >■ oo limit. This is because in the limit g'' — )■ oo, the 
product of operators appearing in two-point correlators can be expanded, via the OPE, in a set of local operators 
with coefficients that are calculable in the few-body sector, in which the theory is integrable, without reference to 
the many-body ground state. While this limit is not directly related to hydrodynamic response at large length and 
time scales, it can be connected to hydrodynamics by employing dispersion relations, as was done in ref. jL2|. This 
gives mo del- independent predictions in the form of sum rules. In this section we briefly summarize the procedure for 
obtaining such sum rules from the OPE and refer to [1^ for details. 

Starting from the Kramers-Kronig dispersion relations for a causal Green's functions x(a;) (such as the functions 
XL,T{i^,q) introduced the previous section). 



Re xi^) = Pr / 

"'0 



oo ^, ,'2 



dio'^ Im x(w') 



. , f^ duj' Re x(w') , ^ 

ImxH - -2c.Pr / ^^, (23) 



where Pr is the Cauchy principal part, one obtains upon taking the formal limit oj -^ oo 

1 [°° duj~ 

n=0 



Re x{^^^)^-Y. J^;^ / ^^''"Ini X(^')- (24) 



Therefore, frequency averages of the dissipative function Imx(i^) are determined by the high-frequency asymptotic 
behavior of the Green's function. In field theory, x(w) generically has power-law ultraviolet behavior and its is not 
valid to expand the dispersion relations. This can be circumvented |22| by employing a Borel resummation of the 
dispersion relation, which gives 

'^e-''-h^x{^)^Y.^M)-^ (25) 



^0 JO 



r(a) 



where the expansion coefficients aa are fixed by the asymptotic behavior of the Green's function at large imaginary 
frequencies 

x(iw->oo)-^a„(w2)-". (26) 

a 

Here we have used the property that x(*'^) is a real function for real lu > Q |24| . If the free parameter wq is taken 
sufficiently large relative to the scales characterizing the ground state, the sum on the RHS of Eq. (|25|) can be truncated 
after the first few terms. It becomes a phenomenological/experimental problem to determine an optimal value of wq 
for which the integral on the LHS of Eq. (|25p is dominated by the low energy data while retaining convergence on the 
RHS. Note that in general, the series is asymptotic although for conformal theories (such as the theory of fermions 
at unitarity) the radius of convergence is finite. 

At the linear response level, the function x('^) can be expressed as a two-point correlation function of operators. 
Thus the coefficients appearing in the expansion of Eq. (j25l) can be understood using the operator product expansion 
(OPE) pi. I1CHL3 [. For a non-relativistic QFT, the OPE is the statement that a product of local operators has the 
following expansion in the limit of zero separation, 

A{x)B{Q) ^^^ Y. \Sf--''^-''^ U (^) Oa(0). (27) 

Eq. (j27l) is an operator statement which holds inside arbitrary matrix elements. The sum on the right hand side is over 
operators with scaling dimension Aq,, and the Wilson coefficients /« ( |a;| /tj are calculable by taking few- particle 
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FIG. 1: Diagrams contributing to the matrix element of the operator M*^''°'(g) between one-particle states. The symbol 
denotes an insertion of the operator r'-' . 



matrix elements on both sides of the OPE relation^. 

Inserting Eq. (j27p inside a time-ordered Green's functions gives the q^ — ^ oo asymptotic behavior 



G^^iq) = * f d'^+^xe"^-^TA{x)B{0)) - ^ ■ 



1 



;[(d+2)+A„-A^-AB]/2^" \2mUJ 



<t 



(040)) 



(28) 



In the w — )■ 00, fixed (f limit considered in this paper, the Wilson coefBcients c^iz) are analytic functions away from 
multi-particle thresholds. The Feynman Green's function G"^^{q) agrees with the causal correlator away from real 
frequencies, and therefore the OPE reproduces the expansion in Eq. (p6)) . relating the coefficients Oq to the condensates 
(Oa)- These condensates are functions of thermodynamic variables that characterize the many-body ground state 
and thus scale as powers (Oa) ^ 11^"'^ of the Fermi energy. The OPE is valid for cj ^ /i (but smaller than the energy 
scale associated with finite effective range corrections to Eq. ([IJ). 



III. VISCOSITY SUM RULES 



A. OPE 



We now determine the asymptotic w — >■ 00 behavior of the frequency-dependent viscosities 77(w), C('^)- In order to 
do so, we construct the OPE for the product 



M'^ 



M 



{u,q = Q)^Jd''+' 



xe'«-"T{r*^(2;)T'^'(0)}, 



(29) 



where t^ (x) is the traceless part of the stress tensor in Eq. (1221) . t^ = T^ — SijTkk/d. For our purposes it is sufficient 
to take g'^ = (w — >■ 00, q^ 0). 

The OPE is constructed by taking matrix elements of Eq. (1^^]) between few-particle states. In the limit q^ ^ oo 
these matrix elements can be expressed by local operator matrix elements, with q dependent coefficients that we 
calculate below. At leading order in the OPE expansion, the necessary matrix elements are between one-particle 
and two-particle states. We can take the limit in which these external states have vanishing energy and momentum. 
This omits the contributions from local operators constructed from at least one pair of Fermi fields {^{x),il'' (x)) 
and at least two derivatives (by rotational invariance, operators with one derivative will ultimately give vanishing 
contributions to Green's functions in the many-body ground state). 

The one-particle matrix element of M'^^''^^{q) is given by the graphs in Fig. [TJ These are exact to all orders in g. 
Given that 



{p'\n,{q)\p) = {2ttY+H{p 



1 

2m 



Pt{p-Q}] +Pj{P~ Q}t 



SijP- (p-q) 



(30) 



^ We are simplifying the discussion slightly by ignoring possible renormalization group scale dependence of the Wilson coefficients (in 
even spatial dimensions), or dependence on scattering length away from the fixed point a — > 00. However, our results in the next section 
take account of the exact a-dependence of the OPE. 






FIG. 2: The exact 2 — >■ 2 scattering amplitude. 






FIG. 3: All tree graph topologies contributing to the the 2 — s- 2 matrix element of the operator M'^''''{q). The shaded blobs 
denote exact vertices, given in Fig. [2] The symbol (g) denotes an insertion of the operator r'-' . 



we see that both diagrams are zero in the limit of zero external momentum. Consequently, there is no contribution 
from operators constructed from one pair {il!{x),ip^ (x)) to the asymptotic behavior of the frequency dependent shear 
viscosity. 

There are non-zero contributions from the two-particle matrix element of A'P^'^^{q). We shall compute those to all 
orders in the coupling g. In order to do so, we require the exact 2 — >■ 2' scattering amplitude, which is shown in Fig. [2] 
and given in d spatial dimensions by 



A-\q) 



(47r)''/- 



,r(l - d/2) 



-m(w- -Eij- 



rf/2-l 



(31) 



where q — (w, q) is the total energy and momentum of the initial state. In particular, for d — i this gives the relation 



mg_ 

An 



(32) 



between scattering length and the coupling constant g, indicating that in dimensional regularization, the unitary limit 
a — >■ oo corresponds to g — >■ oo. 

The full set of diagrams contributing to the 2 — > 2' matrix element are given in Fig. [3] and Fig. |31 In the limit of 
zero external momentum and g — ?> 0, the diagrams with operator insertions on external lines (all graphs in Fig.[3l and 
graphs (a) in Fig. U) vanish identically. It follows from Eq. ([50]) that each loop integrand in Fig. HJd) is proportional 



to edj 



Sij (where £ is loop momentum). By rotational invariance, the angular integrals give IJj — >■ Sijl /d and 



d"'J. 



therefore Fig. IH^d) vanishes. 



The remaining diagrams are, after performing the energy integrals by residues, 



2 r ^df 1 1 I Tfi 



Fig.sib) - ^/ (^^.^^r^^^l ^^^^^. - "-M V^'^ - "-M + (- - --)' (^^) 



and 



iv?u: J (2^)'i R-w - 2E.+ le \ ^'^' d ^'' ] ( ^^'^^ rf 



Fig. He) = ^^ / j,-^- ,, ,^ , ,^ ( ^./, - -T<5., ) ( idi - -Ski 1 + (c^ ^ -c^), (34) 



n 



where Er— I /2m. Using rotational invariance, we make the replacements 



d' 



Idj -^ -jSij (35) 



and use 



d^£ E^} 



kljlkk -> ,. , , [5tj5ki + 5rk5ji + 5^l5kj] , (36) 



1 M ^" ^ r(l-n-rf/2)r(d/2 + n) ^ . ^ . .^d/s+n-i .„„. 



This gives the exact 2 — > 2' matrix element at zero external momentum and q ~ Q: 

(2'|M->'(g)|2) = -^i^^-^J^^ [(-™- - ^^f'-' + (mu - ^ef^-'] ■ (^,,5,, + 5,4, - ^<J.,4.) , (38) 

which indicates that the OPE for M^^'''''{q — ^ oo) contains the four-fermion contact operator O4 = '0{'02V'iV'2- Using 
(2'|C'4(0)|2) ^ ^A{0)^/g^ (see 0), we obtain 

M^«-n- ^ 00, <f = 0) ^ ^(^7^) ^ll^O^y ^ [^""" " "^'^'"' ^ ^"" " "^'^'"'1 ■ ^^'"^'^^^^ ^ " ■ ' ^^^^ 

where . . . denotes the contribution to the OPE from operators of higher dimension. 

B. Asymptotics and Sum Rules 

The results of the previous section can be combined with those of ref. [12| on the dynamical structure function 
S{uj,q) to obtain the asymptotic behavior of the frequency dependent viscosities. We will focus on the case d = 3. 
By Eq. ([H]), and Eq. 1^ we obtain 

Rertiiu) - ^(TOu;)-i/2.c(a)+of^y . (40) 

where we have used the relation Imx^'^ '"^^ {i^jJlIr = tanh(a;/2T) • Re {M''y'^y{uj,q)) valid for real u > 0. In Eq. pO]) . 
C{a) is the by now ubiquitous Tan contact [l|-|3|, which is the expectation value C{a) = (—g^On) and is proportional 
to the thermodynamic quantity dT/da~^ [20|. 

The fact that Re77(w) ^ u;~^/^ at large frequencies is a simple consequence of Eq. \2d>\ . given that the scaling 
dimensions are A(Tij) = 5 and /S.{g^04) — 4. The numerical coefficient is in agreement with the result of [19[. However, 
it follows from our analysis that Eq. (1401) is exact to all orders in the coupling constant g, and consequently gives 
the exact dependence on the parameter a^/rmj^ for arbitrary values of the scattering length o. There are corrections 
to Eq. (|40|) from operators starting at dimension A = 6 in the OPE (for bosons, there would be corrections from 
three-body contact operators at A = 5 [2Q|). Such terms are suppressed by a factor of (/^/w)^ in the large frequency 
limit. 

Using Eq. (fTT]) . Eq. pO)) and Eq. (JA25P from appendix |X] (reproduced from ref. [l2]) we also get the asymptotic 







FIG. 4; All loop graph topologies contributing to the the 2 — )■ 2 matrix element of the operator M^-''^\q). The shaded blobs 
denote exact vertices, given in Fig. [5] The symbol ® denotes an insertion of the operator r'-* . 



bulk viscosity, 



ReC(t^) 






(41) 



which vanishes in the limit a — > cxd, as expected by Schrodinger invariance [9|. 
From Eq. (ITi|) . and Eq. ([M)) one gets the following sum rule for Re 77(0;): 



■^0 "'0 



00 A, .2 



^'^ e-'^-ouRe^{Lo) 



60tt 



rr(l/4)(mc^o)'/'-C(a)+0 



(42) 



which exhibits the exact functional dependence on the parameter l/ay/rmJo at this order in the OPE. From Eq. (fT7|) 
and Eq. (IA27[) we also obtain a bulk viscosity sum rule 



1 r^ dw2 _„,2,.,,2 „ ^_ V2m 



^0 "'0 



-e-" /"oa;ReC(w) 



{mLJoY^^Fia^T^^) ■C{a)+0 
don \uj 

V2r, 



^^V 



367rr(l/4) 

V2 



{mujo) {a^mcjo) x 



1 - ^r(l/4)2 (aV^^i^) ' - %/2r(l/4) (a^THZJ^) ' + 



■Cia) + 0[^y 



(43) 



where the function F(a^mujQ) is obtained by Borel resumming the Wilson coefficient in Eq. (jA27|) . It is given by 

F(aV7i^) = z'^/^^giz) , (44) 

where 



7(z) = V2(l-e-^)+2 



5/4 



.2%/2, 



j;^^e-^7(5/4, -z) - zV4^r(l/4)e-^7(7/4, -z), 



(45) 



and 7(s, z) is defined by the analytic continuation of the integral 

j{s,z)^ z^' [ dtt'-'^e''^ (Rez>0). (46) 

Jo 

The function 7(5, z) is proportional to e^^ times the confluent hypergeometric function ii^i(l, s, z) and is therefore an 
entire function over the complex z-plane for fixed s. Both sum rules receive corrections from operators at dimension 
A = 6. 

IV. CONCLUSION 

In this paper we have obtained new Borel sum rules for the spectral shear and bulk viscosity of a non-relativistic 
Fermi gas. By adjusting the Borel parameter wq, the sum rules can be used to constrain the long time scale transport 
behavior of these systems. Our results are exact in a, but receive corrections of relative size (/i/cjo)^ ^ 1 from 
condensates of dimension A > 5. 

We find it particularly striking that the scattering length dependence of our shear viscosity sum rule is given exactly 
by Tan's function C{a). Thus, at least for the particular case of the shear viscosity, our results provide a simple direct 
link between transport and the theoretically more tractable thermodynamic behavior of this strongly coupled system. 
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Appendix A: Current-Current OPE 

In this appendix we give results for the operator product expansion of 

M'^'{q)= f d'^+^xe"^''T{r{x).P{0)}, (Al) 

in the limit oj — ^ 00 but arbitrary (^ j^raio. Our results are valid up to corrections from operators of dimension A > 5. 
For terms up to operator dimension A = 5 we find exact Wilson coefficients to all orders in the coupling constant g. 
We take as our starting point the OPE for the density correlator, 

M{q) = I d'^+ix e^9-^ T {n(x)n(0)} , (A2) 

which was computed including the contributions of operators A < 5 in ref. [12|. As we use that result in the 
determination of the AP'^{q) OPE, we now recall them. Matching in the one-particle sector generates the operators: 

n, din, Ji, didjU, {diJj +djJi), 



^ v^t^^^, V't ^,d^ + Y_y^ ^^ (^_,a. + l-y, 



(A3) 



in the OPE of M{q), while matching in the two-particle sector gives rise to the additional operator 

Oi = ^14^1^2. (A4) 

The expression for the OPE simplifies if one restricts attention to matrix elements in a translationally and rotationally 
invariant many-body state \Q). In that case, excluding n, the VEVs of the operators appearing in the first row of 



10 

Ea. (IA3l) vanish. The resuh is then ((•••) denotes a matrix element in the state \il), or a thermodynamic average at 
finite temperature), 



{M{q)) = Cniq) (n) + Cu{q) (H) + C,{q) (g^O^) 



(A5) 



where 



H = — \VM^ - g44i;i^2, 



(A6) 



Cn{q) = 



[lo - E^){uj + E^) 



(A7) 



CH{q) - -E^ 



7 



(uj - S,-)3 {uj + E^f 



(A8) 



C,{q) -* 



A{q)(h{q)--— ^' ''" ' 



.g a; - i^,- 



T^hiq) - -h{q) + {q ^ -q) 



2i 
9 



oj — E^ Lo + E^ 



4dE^ 
9 d 



1 



{lu - E,~r (^ + E^r 



(A9) 



and 



Uq) 



I r(a + i-f)r(f-a) 



V47r/ 



r(i) 



^-Eq' — uj — ie 



. d d 

2F1 [a,a + l- -,-, 



5^? 



2' 2' uj~^E^+ie 



(AlO) 



Note that 



/o(g) = - ^ 



5 ^(9)' 



(All) 



where A{q) is the full 2-^2 scattering amplitude from Eq. (|3ip , and that the integral in Eq. (|37p corresponds to the 
case I-n{u},q = 0). 

In order to compute M^'^{q), we have found it convenient to use Ward identities (current conservation) to relate it 
to the OPE of M{q) given above and to the OPE with one insertion of the current J and one insertion of the density 
operator n. We only display results for the expectation value (Af*'^(g)) in a state that is invariant under rotations 
and translation, as the full operator relation is not useful for obtaining the physical observables discussed in sec. HUB] 
In this case, we find 



{A,P^^iq))^l^L{q)+[S^^- 



q'qJ 



Tiq), 



(A12) 



where 



L{q) = Ci'\q) (n) + &^\q) {%) + Cl}\q) (g'O,) 
T{q)^d^Hq){H)+CPiq){g'0,), 



(A13) 
(A14) 
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and 



CW(g)^^C„(g)--, 
g^ m 



(P 






E? 



d m{uj — Eq){u) + E^) ■ 



(A15) 
(A16) 
(A17) 
(A18) 



C?\q) 



2m 



{d - l)ujEij 



I-M + 



3to 



[d - 1)£;,- 



h{q) - 



(^.- 



m{Ei- 2ujEff + 3w^ 



m(uj - Eff) ^ , , , , 



^lE^ 



2{d- 1)ujE^ 
1 



Mq) 



g d m(uj ^ Eq){uj + E^) 



The causal and time-ordered Green's functions for the current were defined in Eq. Q and Eq. 
Since G^'^{q) — iAP'^{q) and IniG'*'^(q) = Imx^R (q) (for real a; > 0), we have for q^ ^- cx) 

Im XL{q)= Re L{q) = {g^Oi)ReCi''\q) = ^ {g^Oi)ReCA,q)- 

In the limit g — > the relevant Wilson coefficients have the expansion, for w > 0, 

4 



ReCi'\q) 



<f 



(mw)2 d{d + 2) 



Im [Aiio)-^ + Ai-u)-'] + O {(f) , 



(A19) 
respectively. 

(A20) 
(A21) 

(A22) 



ReCciq) 



^ 



mu) J d'^{d + 2)uj'- 



-Im 



AH 



8[d-l)-{d-2f{d + 2) 



A{u) 



+ (w -> -uj) 



+ O {<f) , 



(A23) 



where A{uj) = A{uj,q — 0) is the 2-^2' amplitude of Eq. pip at center-of-mass energy ui. In d = 3 spatial dimensions, 
these expressions become 



m q 



Re C(2)(,) = -^{^] {mu)-^n + q {t) 
loTT V mui ' 



(A24) 



and 



Re CM 



Am ( q 



nnu!) 



-3/2 



16 (1 + (aVmZJ)2) 



+ o{f). 



457r \muj 
The Borel sum rules are related to the Wilson coefficients at imaginary frequency. In d = 3 we obtain for ujq > 0, 



(A25) 



lOTT \mujQ J 



5^/2 1 

1 —[ay'mujo) 



o{<f) 



(A26) 



Cc(za;o) 



7V2im^ f q^ 



GOtT \77lWo 



(mwo)' 



-3/2 



21 (l - V2(a^/?H^)-i -F (a7?HZJ^)-2) 



O {cf) , (A27) 
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